Self-reptation and slow topological time scale of knotted polymers 
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We investigate the Rouse dynamics of a flexible ring polymer with a prime knot. Within a Monte 
Carlo approach, we locate the knot, follow its diffusion, and observe the fluctuations of its length. 
We characterise a topological time scale, and show that it is related to a self-reptation of the knotted 
region. The associated dynamical exponent, zt ~ 2.32 ± .1, can be related to that of the equilibrium 
knot length distribution and determines the behaviour of several dynamical quantities. 
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In the physics of polymers, mutual and self entangle- 
ments play an essential role [l|]. A particularly relevant 
type of entanglement is associated with the presence of 
knots. Indeed, it has been known for some time that 
long ring polymers inevitably contain a knot 0. Topol- 
ogy is also of much interest for biopolymers, where knots 
have been found in the DNA of viruses and bacteria 
and also in some proteins Q- The recent experimental 
possibility to tie knots in DNA double strands or actin 
filaments , and to study their diffusion in the presence 
of a stretching force @ , has further increased the interest 
in topology related issues among polymer physicists. 

So far, the statistical physics of knotted polymers has 
focused mostly on static properties 7]. For example, it 
was found that the presence of a knot does not alter the 
exponent ly that relates the radius of gyration Rg{N) to 
the length N of the polymer, i?p(iV) ~ iV, but only in- 
fluences scaling corrections ^ lOj . There have been fewer 
studies of dynamical scaling properties of knotted poly- 
mers, even though the associated time scales could be rel- 
evant to describe gel electrophoresis, folding of knotted 
proteins or other dynamical processes. Simulations have 
given evidence that a peculiar, topological, characteristic 
time tt determines the decay of the radius of gyration 
autocorrelation function of such polymers l3, Hi- This 
decay time appears longer than that observed for open or 
closed unknotted chains. As a function of chain length, 
it was found to scale with a dynamical exponent zt, that 
could not be determined precisely but whose value prob- 
ably is bigger than the Rouse one. So far, the physical 
origin of this time was not understood. 

In this Letter, we investigate how a knot influences dy- 
namical scaling properties of flexible polymers in a good 
solvent within the framework of Rouse dynamics [H Il2l| . 
In a simulation, we directly follow the motion of the knot- 
ted region and the fluctuations of its length (measured 
along the backbone). We develop a simple picture, based 
on reptation, that establishes a link between zt and the 
exponent governing the distribution of this length. 



In Rouse dynamics, one models a polymer as a set of 
A^-|-l beads (monomers) that are connected by harmonic 
springs, and that are subjected to random thermal forces 
exerted by the solvent. The motion of the monomers is 
described by a Langevin equation, which for ideal chains 
can be solved using a transformation to normal coor- 
dinates. If self-avoidance is taken into account, Rouse 
dynamics can no longer be solved exactly, but scaling ar- 
guments together with numerical results [l^l, show that 
the center of mass of the polymer, i?cm, performs ordi- 
nary diffusion, i.e. 

g3it,N) = {(^Rcmit) - i?c,„(0))') ~ 6D„r^t (1) 

where the average is taken over realisations of the 
stochastic process. The diffusion constant is inversely 
proportional to the length N of the polymer, i.e. Dcm ~ 
N^^. The autocorrelation function of the end-to-end vec- 
tor decays exponentially with a time constant tji that 
grows as 
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with N. In d 3, 2i/ + 1 ~ 2.2 [13|. We will refer 
to tr as the Rouse time scale. It can be interpreted as 
the time that the center of mass of the polymer needs 
to diffuse over a length equal to its radius of gyration, 

Tfl~i?^(iV)/D,„~iV2.+ l^ 

The diffusion of one particular monomer, or of a seg- 
ment of TO + 1 monomers, is described by the function 
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It is known 14] that gi{t) has a scaling form 
g,{t,N)^tPF{t/Tn). 



(4) 



The function F describes the crossover between an initial 
power law regime {F(x) — > constant, for x — > 0) and a 
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late time regime for which a monomer has to follow the 
diffusion of the whole polymer as given by ([1]). This 
implies a power law form for F at late times and P — 
2^/(2j/ + 1), i.e. in the initial regime the movement of 
a segment of the polymer is subdiffusive. Notice that 
gi{TR,N) ~ Rq{N), so that the Rouse time scale also 
corresponds to the time for one monomer to diffuse over 
a distance equal to Ra{N). In the whole dynamics, it is 
the only relevant time. 

While Rouse dynamics neglects important physical ef- 
fects, such as hydrodynamic interactions between the 
monomers, several experimental situations are known by 
now for which it gives an adequate description. As an ex- 
ample, we mention that current fluorescence techniques 
allow to follow the motion of individual monomers in, 
e.g., DNA-chains and hence to directly determine a func- 
tion such as gi{t, N). Measurements of this type on dou- 
ble stranded DNA have recently seen diffusive behaviour 
of individual monomers consistent with Rouse behaviour 

A crucial quantity to characterise the dynamics of a 
knotted polymer turns out to be the length of the knot. 
A precise definition for this quantity can be given for flat 
knots [3l ■ These are knots in a polymer that is strongly 
adsorbed to a plane or constrained between two walls. In 
this context, it was found that the length Ik of a knot k 
is a fluctuating quantity, whose equilibrium distribution 
function is a power law, p{lk) ^ ^k'^G{lk/N) where G is a 
scaling function. It follows, that the average length of the 
knot scales with N as {Ik) ~ -/V^, with a = max[0, 2 — c]. 
In a good solvent, flat knots were found to be strongly 
localised (ct = 0) [16i|, while below the 0-transition, they 
delocalise [a = 1) [17]. In order to extend these results 
to three dimensions, one needs a good definition of the 
length of a knot. This should correspond to the intuitive 
idea that it measures the portion of the poly mer back- 
bone which "hosts" the knot entanglement [19| . Recently, 
a powerful computational approach to determine such a 
length, and its scahng properties, was introduced In 
this method, for a given knotted ring polymer, various 
open portions are considered and for each of these a clo- 
sure is made by joining its ends with an off-lattice path. 
The length of the knot in a given configuration can then 
be identified with the shortest portion still displaying the 
original knot. In this way, it was found that in good sol- 
vent, three dimensional knots are weakly localised with 
an exponent tr — .75 ± .05. 

To simulate the dynamics of knotted polymers, we 
start from a simple self-avoiding polygon (SAP) on a cu- 
bic lattice with a knot in it. Most of our calculations 
were done with a trefoil (3i) knot. After performing a 
number of BFACF [2^ moves to relax the configuration, 
the resulting SAP evolves according to a iV-conserving 
Monte-Carlo dynamics with local moves only. This is ex- 
pected to give a proper description of Rouse dynamics 
[2l|. During the subsequent evolution, using the compu- 



tational techniques developed in [9|, observables such as 
the length of the knot, the radius of gyration of the whole 
polymer, the location of the center of mass of the poly- 
gon and the location of the center of mass of the knot, 
Rcm,k, are computed. We have performed calculations 
for various N < 400. 

From the data on the knot length as a function of time, 
we can construct the knot length autocorrelation func- 
tion {lk{t)lk{0))c (where the subscript 'c' indicates the 
properly normalised connected autocorrelation). In Fig. 
[U the top (red) line shows our results for this quantity 
for N = 400. This autocorrelation function has a simple 
exponential decay. Similar behaviour is found for other 
N values. Fig. [2] presents our data for the decay time 
constants as a function of N. The behaviour is power 
law, and the exponent is 2.33 ± 0.08. The value of this 
exponent is higher than that of Rouse dynamics. Further 
evidence that this presents a new time scale comes from 
an investigation of the autocorrelation function of the 
radius of gyration of the polymer. In Fig. [T]we also plot- 
ted our results for this quantity (bottom hnc). Clearly, 
it has a double exponential decay. In comparison, for an 
unknot, we observe a pure exponential decay. A care- 
ful analysis shows that the time constant of the unknot, 
together with the fast one of the knotted polymer, are 
proportional to tj^. On the other hand, as can be seen 
in Fig. [l] the late time decay of the radius of gyration 
autocorrelation occurs with the same time constant as 
that of the length autocorrelation. A quantitative analy- 
sis of its A^-dependence supports this conclusion (see Fig. 
[5]). Indeed, the value of the associated exponent equals 
2.31±0.08, consistent within the numerical accuracy with 
that determined from the length autocorrelation. To- 
gether, these data therefore provide strong evidence for 
the presence of a new, slow time scale tt ~ N^^, with 
zt = 2.32±.l, in the dynamics of knotted polymers. The 
data on the knot length autocorrelation show that this 
scale corresponds to the time over which the knot length 
decorrelates. 

To check if the results are influenced by the topology of 
the knot, we have performed a completely similar study 
for the figure-eight (4i) knot. Again the radius of gy- 
ration autocorrelation function decays as a double ex- 
ponential, with a fast time scale that can be identified 
with Tji. The slower time scale grows with N with an 
exponent whose numerical value, 2.33 ± .08 is consistent 
with that found for the trefoil. Moreover, the same ex- 
ponent was found to govern the decay of the 4i length 
autocorrelation function. 

A simple argument explaining the value of zt is based 
on a reptation picture. Reptation describes the motion 
of a polymer in a dense environment of other polymers 
that restrict its movement to a tube P, This is, 

for example, the case in a polymer melt. Through a 
diffusive motion along the tube axis, first one end and 
later the whole polymer moves out of the original tube, 
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FIG. 1: (Color online) Semilogarithmic plot of the autocorre- 
lation functions versus time for a SAP with = 400. Shown 
are results for the radius of gyration (black, below) and the 
knot length (red, above). The blue line was obtained by shift- 
ing ' ' 



(5) 



This function has some similarity to gi{t^ N) since it de- 
scribes the motion of a segment of the polymer. Hence 
we can expect it to have a scaling behaviour similar to 
There are however differences since the number of 
monomers in the knot increases with iV^ on average, and, 
at fixed N , fluctuates in time. 

In Fig. [31 we show a typical result {N = 250) for 
the functions g2{t,N) and g3,k{t,N). Our data for the 
diffusion of the center of mass of the whole polymer are 
fully consistent with Eq. Il} and with the relation Dcm ~ 
N^^, a strong indication that the diffusion of the polymer 
as a whole is not changed by the presence of the knot 11 1 . 
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FIG. 2: (Color online) Log-log plot of the topological time 
scale versus N as determined from the radius of gyration au- 
tocorrelation function (black circles) and from the knot length 
autocorrelation (red squares). The fitted lines have slopes 
-2.31 ± 0.08 (black dashed line) and -2.33 ± 0.08 (blue full 
line) respectively. 

creating in this way a new one. This process is slow and 
it takes a time ~ for the new and the old tube to 
decorrelate. In the case of interest here, we can imagine 
that the knot, which is weakly localized, creates a local 
tube that constrains the movement of the monomers it 
contains. Since it is only the knotted part, of length 
proportional to N" , that performs the reptation in this 
case, we expect a decorrelation of the knot tube after 
a time N^"'. Taking the estimate of a from [S], we find 
3(T ~ 2.25±.15, which is consistent with the identification 
zt — 3cr. This agreement supports the plausibility of the 
self-reptation mechanism. 

In order to determine whether the topological time also 
influences other dynamical properties, we have investi- 
gated the diffusion of the center of mass of the whole 
polymer and of the knotted region. From our data, we 
determine (73 (t, N) and a similar quantity for the knotted 
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FIG. 3: (Color online) The function of the whole polymer 
(lower curve) and of the knot (upper curve) (A = 250). The 
straight lines are best fits through the initial (blue) and inter- 
mediate (red) time regime of the knot. Their slopes are .29 
and .78, respectively. The dashed vertical line correspond to 

More interesting is the diffusion of the knot itself. As 
can be seen in Fig. [31 gj,^^ shows several distinct power 
law regimes. The value of the estimated exponent for the 
early time region fluctuates with N , but a good average is 
0.27 ±0.1. This initial regime ends after a time ri, which 
grows as a power of A^, t\ ~ A^^i, where z\ = 1.97 ± .1. 
After a short crossover, in an intermediate time regime, 
the behaviour is again power law, 173^^ ~ t'' . The value of 
7 decreases from « .81 for A^ = 100 to a value close to .66 
at A^ = 400. An extrapolation gives the asymptotic value 
7 — .6 ± .03. Finally, and in analogy to the behaviour 
of 51, we expect a crossover of gs^k to linear behaviour, 
since the knot eventually has to follow the whole polymer. 
This crossover has not happened yet for the times we 
were able to simulate. In Fig. [31 the vertical dashed 
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line indicates the Rouse time, which we have estimated 
from the relation g3{T]i,N) = Rq{N). We checked that 
the Rouse time calculated in this way indeed grows as 
N"^-"^. We thus see that, as was the case for the length 
autocorrelation, tr doesn't play a special role for g^^k, 
and moreover, it seems that the second crossover occurs 
on a much slower time scale, which we expect to be tt- 
In analogy with ([4]), we therefore propose this second 
crossover to be of the form 



g3,k{t,N)^t-<H{t/TT) 



t > Tl 



(6) 



where H{x) becomes constant for small x. Since for t > 
Tt, the behaviour of ([6]) must cross over into that of ([T]), 
we obtain the relation (1 — j)zt — 1- Using the estimate 
zt = 2.32 ± .1, this leads to 7 = .57 ± .02 consistent 
with our estimate. In Fig. [H we show a scaling plot 
of our data for g3,k{t,N) for various iV- values (leaving 
out the initial power law regime) and using the values 
Zt = 2.35,7 = -57. The scaling is rather well satisfied. 
So, the available numerical evidence is consistent with 
the identification of the second crossover time with tt- 
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FIG. 4: (Color online) Scaling plot of our results for 
g3,k/N'''^'' versus t/N^'^ for various A^'-values using zt ~ 
2.35, 7 — .57. Different colors correspond to different A'^ val- 
ues, red: N =150, blue: 200, dark green: 250, yellow: 300 
and light green: 350. 

In conclusion, by the first calculation in which dynami- 
cal information on the knot itself was monitored, we pro- 
vided convincing evidence, that the presence of a knot 
in a ring polymer introduces a new slow topological time 
scale that is due to a self-reptation of the knotted re- 
gion. Our results thus show the physical importance of 
the fact that the length of a knot in good solvent is a 



scale invariant, fluctuating quantity. 

The picture presented here implies that for a polymer 
below the 0-point, where the knot delocalises (cr = 1) 0], 
the topological time becomes proportional to the repta- 
tion time of the whole globule N^. Moreover, in 

cases of strong localization, as for flat knots, we expect 
that the extra, slow time scale should not be present. 
The verification of these predictions and the experimental 
search for the topological time scale in specific processes 
are challenges for future work. 
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